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ABSTRACT 

We present a theory for interpreting the sodium lines detected in transmission spectra of exoplanetary at¬ 
mospheres. Previous analyses employed the isothermal approximation and dealt only with the transit radius. 

By recognising the absorption depth and the transit radius as being independent observables, we develop a 
theory for jointly interpreting both quantities, which allows us to infer the temperatures and number densities 
associated with the sodium lines. We are able to treat a non-isothermal situation with a constant temperature 
gradient. Our novel diagnostics take the form of simple-to-use algebraic formulae and require measurements of 
the transit radii (and their corresponding absorption depths) at line center and in the line wing for both sodium 
lines. We apply our diagnostics to the HARPS data of HD 189733b, conhrm the upper atmospheric heating 
reported bv iHuitson et al.l (120121) . derive a temperature gradient of 0.4376 ± 0.0154 K km“^ and hnd densities 
^ 1-10"^ cm“^. 

Subject headings: radiative transfer — planets and satellites; atmospheres — methods; analytical 


1. INTRODUCTION 

Due to its large cross section and favorable wave¬ 
length range, the sodium doublet lines have been a boon 
to astrono mers seeking to characte ri se exoplaneta ry at¬ 
mospheres (ISeager & Sasselo^ l2000t iBrow 3 UMl- In 
fact, the hrst detection of an exoplanetary atmosphere 
(that of HD 209458b) was accomplished via measur¬ 
ing the sodium doublet (ICharbonneau et akl l2002h . Ever 
since that discovery, sodium has been detected in several 
hot Jupiters using both space- and ground-based transmi s- 


sion spectroscoDV ( 

Redfield et al.l 120081; iSin? et alJl2008ayhl: 

Snellen et al.l 12008 

; Vidal-Madiar et all 1201 U iWood et a1 

2011 

Huitson et al.l 120121; Sins et alJ 20121 Zhou & Bavliss 

2 OI 2 I 

PontetalJ 2013 

; iNikolov et al.l 120141 iBurton et al. 

2 OI 5 I 

iWvttenbach et al.l 

2015|1. 


As better instruments come online and our ability to resolve 
the sodium lines improves, it is worth revisiting and redevel¬ 
oping a theory of how to interpret them. From a remote sens¬ 
ing perspective, a pair of fully resolved sodium lines can, in 
principle, map out the temperature-pressure prohle of an at¬ 
mosphere at high altitudes and yield the sodium abundance. 
A qua ntity used to interpret absorp tion lines is the equivalent 
width (ISpitzejl 1 9781 ; lDrainell20 1 lli . 

W = [ Wx dX, 

J ( 1 ) 

Wa = 1 - exp(-r), 

where t is the optical depth of the intervening material and 
A denotes the wavelength. The equivalent width is essentially 
the width of a box with the same depth as the trough of the 
absorption line, such that it encompasses the same area. 

While it is not immediately obvious, Wx and the transit 
radius (R) are actually independent observables. A simple 
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Fig. 1.— Schematic of transit radii associated with the line center and wings 
of the sodium doublet. As each line is expected to have the largest cross 
section at line center, the corresponding transit radius (i?o) is the largest. In 
the line wings, where the cross section is the smallest, the transit radius (H) is 
correspond ingly smaller. The cartoon on the right reproduces the schematic 
depicted in lFortne^ 420051) . 

thought experiment demonstrates this. Consider a hctitious 
star that emits achromatically. At the exact moment of tran¬ 
sit, one may record an absorption spectrum of the exoplan¬ 
etary atmosphere. If the wavelength coverage is sufficient, 
the sodium doublet and the continuum and thus Wx may be 
measured. To record the transit radius requires temporal in¬ 
formation: the change in flux in and out of transit. In short, 
Wx measures changes in absorption across wavelength, while 
R derives from the change in flux across time. Certainly, stars 
are not achromatic light sources and one needs to measure the 
stellar lines in and out of transit to properly subtract out their 
influence, but this is an observational, rather than a theoreti¬ 
cal, obstacle. 

Nevertheless, a theoretical challenge with interpreting W, 
for transmission spectra of exoplanets, is that the optical depth 
depends on the transit radius and the radius itself depends on 
wavelength (Figure [T]). It is generally difficult to derive the 
functional form of i?(A), but it is straightforward to measure 
it. Additionally, while W is associated with a fixed sightline 
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in traditional studies of the interstellar medium, it is associ¬ 
ated with a set of sightlines in transmission spectra, each cor¬ 
responding to a different wavelength. Instead of a theory for 
W, we develop one for W\{R). Observationally, it takes the 
form. 


VkA = 


F,-F 

Fc 


( 2 ) 


where Fc is the flux associated with the continuum and F is 
the flux at any point within the line. More specifically, both 
Fc and F are integrated over a (small) wavelength interval 6X 
that is chosen based on practical constraints. We shall term 
F/Fc the absorption depth (and leave W\ nameless). 


2. THEORY 


2.1. Order-of-Magnitude Estimates 

Since the line-center wavelength of both sodium lines is 
Ao ~ 0.6 pm, the line-center frequency is vq = c/Aq ~ 
5 X 10^^ Hz. This means that the Doppler line width is 




« 5 X 10® Hz 


(liJit) 


1/2 



( 3 ) 


where c is the speed of light, fee is the Boltzmann constant, T 
is the temperature, m is the mean molecular mass and mn is 
the mass of the hydrogen atom. We assume that sodium is a 
trace element in the atmosphere and that Doppler broadening 
is mediated by a dominant buffer or inert gas (i.e., molecular 
hydrogen). 

By contrast, the natural line width is 

rL = 4^«10^Hz, (4) 

ZTT 

where A 21 ~ 6 x 10^ s“^ is the Einstein A-coefficient. Fur¬ 
thermore, the damping coefficient associated with the Voigt 
profile is 


_ Tl 

2rD 


10 “® 


103 Ky 


-1/2 



( 5 ) 


implying that the Voigt and Doppler profiles are essentially 
identical near line center. 

If pressure broadening is significant, then El needs to be 
replaced by El +Econ in the Lorentz profile, where Ecoii is the 
collisional frequency (lMihalaslll970h . Pressure broadening is 
similar in nature to Doppler broadening. 


coll 


ApEp 

^mfp 


= Perm 


(—f 

\mkBT ) 


(6) 


where /mtp is the mean free path of collisions and (Th 2 
10“35 cm® is the cross section of collisions with hydrogen 
molecules. Pressure broadening may be ignored if the pres¬ 
sure being sensed is approximately 


P < 1 mbar 


m T \ 
2mn 103 K ) 


1/2 


( 7 ) 


We will neglect pressure broadening for our analysis of the 
sodium lines, while being aware of a shortcoming of our the¬ 
ory: we cannot directly diagnose the value of the total pres- 



Fig. 2.— Normalized Voigt, Doppler and Lorentz profiles for ap = 10“3 
and 0.1. The former value is representative of the sodium lines in hot ex¬ 
oplanetary atmospheres, while the latter value is presented as an extreme 
case. The line core is well approximated by a Doppler profile out to sev¬ 
eral Doppler widths. Beyond this, the line wings are better approximated by 
a Lorentz profile. From inspecting these profiles, it is clear that the shape and 
normalization of the Doppler core are insensitive to temperature, whereas the 
difference in normalization between it and the Lorentzian wings serves as a 
temperature diagnostic. 


sure being sensed, since it is degenerate with the sodium abun¬ 
dance. Non-Lorentzian corrections to the sodium line wings 
are only impor tant far away (> l OOEn) from line center and 
for P > 1 bar (lAllard et al.ll2Q^ . 

Consider a point on one of the sodium lines that is separated 
by A A, in wavelength, from line center. In frequency, the 
separation is Az/. If Az/ = Ed, then we have 


AA 


A®Ai/ 

c 


0.06 A 


103 k 


1/2 



( 8 ) 


implying that the lines are tens of Doppler widths wide. Even 
by moving 0.3 A away from line center, one is already well 
into the Lorentzian wings of the profile. In the line wings, we 
may approximate the line shape with a Lorentz profile. 

Collectively, the sodium line profile in hot exoplanetary at¬ 
mospheres (T ^ 1000 K) is well-approximated by a Doppler 
core and Lorentzian wings (Figure |2]l, which allows us to cir¬ 
cumvent the more challenging task of inverting a Voigt pro¬ 
file. 


2.2. Review of Previous Analytical Formulae 

In iFortnevI (120051) . the assumptions of hydrostatic equilib¬ 
rium and an isothermal atmosphere led to an expression for 
the number density. 


n = n-cef exp 



( 9 ) 


where n^ef is a reference value of the number density, z is 
the height above some reference radius R, H = fceT /mg is 
the pressure scale height and g is the surface gravity. Using 
Pythagoras’s theorem and assuming z <C R. iFortneM (l2005h 
obtained (see Figure[T]for the geometry) 


Z Ri 


m' 


( 10 ) 
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It follows that the optical depth of the chord associated with 
the transit radius is 


dx = rtrefcr , 

( 11 ) 

where a is the absorption cross section. In hindsight, one 
could have obtained a quick answer by writing r = n^ei<xX 
and recognising that X ^ v /HR is roughly the geometric 

mean of two vastly different length scales, H and R. _ 

Seve ral properties associated with the formula of iFortne"^ 
(l2005li . restated in our equation (fTITl . are worth emphasizing. 
First, it can only be applied to isothermal situations. Second, 
i? is a wavelength-dependent quantity. Third, since one is 
sensing different values of R across wavelength, the value of 
TT-ref Sampled varies as well. The number density can be fixed 
by writing nj-ef = n-J-gf exp {—z/H), where is its value in 
the line wings. Unlike riref, is a wavelength-independent 
quantity. 

The third line of reasoning allowed 
Lecavelier des Etangs et al.] (|2008|) to invert the expres¬ 
sion for T and obtain 


(*+oo 


T — 


exp - 


2HR, 


z 


iJ In 


T 


( 12 ) 


Since the transit radius always picks out r ^ 1 by definition 
and T is not a wavelength-dependent quantity, it follows that 


mg dz r()(ln(T) lc)(lnii) 

^ [ d\ ^ 2 d\ 


(13) 


The second term within the square brackets is typically 
smaller than the first, i.e., d{\nR)/d\ <C d{h\(j)/d\, since 
we expect the transit radius to vary by fractions of a percent, 
whereas the cross section is expected to vary by orders of 
magn itude. Dropping the second term reproduces equation 
(2) of iLecavelier des Etangs et ^(|2008|) . Equation (fOT l is di¬ 
rectly applicable to situations where dzjdX is constant (e.g., 
Rayleigh scattering), but is less amenable to analyzing spec¬ 
tral line shapes, where dzjdX is non-constant by definition. 


2.3. Obtaining Temperature and Density in the Isothermal 

Limit 

At lin e center, the absorption cross section of the sodium 
atom is (lDrainell2011h 


e^/iuAp / 7rm 
meC \2kB,T) 


(14) 


where /lu is the oscillator strength. The mass and charge of 
the electron are given by mg and e, respectively. 

By denoting the line-center value of W\ as Wxg , we may 
express the line-center optical depth as 

To =-ln(l - VFao) =-In (15) 


where Fg and are the flux at line center and of the contin¬ 
uum, respectively. For the remainder of the study, we shall 
deal only with the absorption depth as an observable. It fol¬ 
lows that the number density associated with line center is 


/ TOgC / _g\ 

\Fj TreViuAo \rJ 


(16) 




Fig. 3.— Sodium line diagnostics adopting m = 2mH, Q = 10^ cm s 
Rq = cm. A/? = 10^ cm, = 0.5 and T' = km~^ for illustra¬ 
tion. Top panel: the line-center number density. Bottom panel: the line-wing 
temperature. For the non-isothermal cases, we set b = 2.5 {dTjdz > 0) 
and 3.5 {dTjdz < 0). 


In the line wings, the optical depth is given by 


r = TT-o exp 



(27ri7i?) 


1/2 7re^/iu$ 
rrieC 


(17) 


with the line shape ($) being approximated by a Lorentz pro¬ 
file sufficiently far away from line center. 


$ = 


^2 iA^Aq 

47rc^ (A — Ao)^ 


(18) 


The transit radius at line center (i?o) in the line wings (i?) 
are related by 

Ro = R + XR. (19) 


Like R, AR is a wavelength-dependent observable. Equation 
(fTTl i can be inverted to solve for the temperature. 


mg AR 
fee In r ’ 


( 20 ) 


by eliminating uq using equation (fThl) . We have defined the 


riQ = — In 
































4 


diagnostic ratio as 

_ In (F/Fc) f irmRo \ AA^ 


■ln(Fo/F,) \2kBRTj 

,3 ln(W) 


An A 


0^21 


i2 X lO-" 


T 


ln(Fo/F,) 
- 1/2 / 


103 K 



r 

\R, 

1 

m ^ 

1/2 


V2mH/ 


AA 

0.1 A 


( 21 ) 


where AA = A — Aq. 

The temperature is mostly controlled by the difference in 
transit radii between the line center and wings. Formally, In r 
has a very weak dependence on T, so equation (l20l l is not an 
explicit formula for T. In practice, one sets T ^ 1000 K, 
measures or assumes all of the other parameter values, uses 
equation (l20l l to compute an updated value of T and iterates. 
A converred answer should obtain within a few iterations. 

Figure|3]shows calculations of no as a function of Fq/Fc for 
both sodium lines. For the sodium Di line, we have used Aq = 
5897.558 A and /lu = 0.320; for the sodium D 2 line, we 
have used Aq = 5891.582 A and /lu = 0.641 (lDrainell201 Ih . 
The Di and D 2 lines are sensing different altitudes within the 
atmosphere. Figure|3]also shows the line-wing temperature as 
a function of the diagnostic ratio. 


2.4. Obtaining Temperature and Density in the 
Non-isothermal Case 

We generalise to a non-isothermal situation by considering 
a linear series expansion of the temperature. 


T = 


dT 

Tref + 

oz 


( 22 ) 


where the non-isothermal pressure scale height is H = 
Tj-ef/T' and the index. 


b = 


mg 

kBT’’ 


(27) 


is the ratio of the non-isothermal to the isothermal scale 
heights. 

For dT/dz < 0, the minimum value of b is set by the adi¬ 
abatic lapse rate such that 6 > (2 -f Noof)/^., where Adof is 
the number of degrees of freedom of the buffer gas. For a 
diatomic gas, we have Noot = 5 and b > 3.5 (ignoring the 
vibrational modes). 

The chord optical depth is 


T = nretCT 





{cos9f^ 


dO, 


where the index is 


Nb 


(2b, 

\ 2 b-l, 


^ > 0 

^ < 0 
dz ^ 


(28) 


(29) 


The integral in equation (l28l l has no general analytical solu¬ 
tion. However, it is soluble if we consider integer values of 
Nb- This is not unre asonable. For example, for HD 189733b 
iHuitson et al.l (1201 2ll estimate T' ft! 1 K km”^^ which yields 
ft 5 for dT/dz > 0. For Nb = 0,2,4 and 6, the integral 
is TT, 7r/2, 37r/8 and 57r/16, respectively. Essentially, the inte¬ 
gral may be represented as ^tt, where C, ft 0.1-1. Numerical 
evaluations of the integral, over a continuous range of values 
of Nb, confirm our intuition (not shown). The optical depth 
thus becomes 


T = C^rrirefcr 




(30) 


where F^ef is a reference temperature and we allow the tem¬ 
perature gradient to be constant and either positive or nega¬ 
tive, 

dT 

— = ±r, (23) 

dz 

where T' is always a positive number. We begin with hydro¬ 
static balance, 

dP 

TT = -P9^ (24) 

dz 

where P is the pressure and p = nm is the mass density. The 
physical assumption being made is that the buffer or inert gas 
(e.g., molecular hydrogen) is in hydrostatic equilibrium and 
the collisional time between itself and the sodium atoms is so 
short that the latter are effectively coupled to the former as a 
single fluid. 

By using the ideal gas law, P = nfeeF, one obtains 



— dn = — 
n 


mg 

fee 


1 


Tre! ± T'Z 


dz. 


(25) 


Here, g = g zL kBT'/m functions like an “effective gravity” 
with an additional term associated with the constant tempera¬ 
ture gradient. Completing the integration yields 


similar in spirit to its isothermal counterpart where we again 
have the characteristic length scale being X ^ s/HR. Gener¬ 
alising from the isothermal situation, we note that Uret, R and 
Fref are wavelength-dependent quantities. 

At line center, we again assume a Doppler profile and derive 
the number density. 


- _ 1 f kBT' \ 

^\FcJ CTre^/iuAo \TrmRoJ 


(31) 


an expression which holds for both dT /dz > 0 and dT/dz < 
0. Unlike before, no now depends also on T'. In the line 
wings, we can again assume a Lorentz profile, eliminate no 
and obtain 


F'Ai? [N/(b+i) - 1 ] \ > 0, 

T'AR[l-r-F(b-i)^-\ ^< 0 . 


The definition for r is identical to the one stated in equation 
(l2Tli . It is worth noting that our expression for F, in equation 
in, does not depend on C- 

Equations (l3Tl l and 1^) constitute an under-determined 
system; two equations and three unknowns (no, F, T'). We 
will break this degeneracy by using measurements from both 
sodium lines to determine T'. 


|r>o. 

[nref(l-0) , |i<0, 


( 26 ) 


3. APPLICATION TO DATA 

To apply the diagnostics to data requires that both sodium 
lines are at least partially resolved, such that for each line we 
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may measure the absorption depths at line center and in the 
line wing. Generally, at the same distance (AA) from line 
center, the line wings of the Di and D 2 sodium lines are sam¬ 
pling different temperatures, which we denote by Ti and T 2 , 
respectively. We denote the line-center temperatures by Tq 1 
and To, 2 - 


1. To determine if the system is isothermal, one computes 
the ratio of the line-wing temperatures. 


T\ ^ ^ 
T 2 Ai?2 ’ 


(33) 


where Ai?i and AR 2 are the differences in transit radii, 
between the line center and wing, for the Di and D 2 
lines, respectively. To a very good approximation, the 
preceding expression holds regardless of whether one 
uses the isothermal or non-isothermal formula for T. If 
Ti = T 2 , the isothermal formulae in equations (fTST l and 
(I 20 I 1 will suffice. 


2. If Ti ^ T 2 , then we first use the isothermal formula in 
equation (l20l l to estimate Ti, T 2 and 


dT _ Ti- T 2 
dz i?i — i ?2 


(34) 


After assessing if the temperature gradient is positive 
or negative, we use the appropriate non-isothermal for¬ 
mula in equation (l32l i to iterate for Ti, T 2 and T'. 


3. With the value of T' in hand, we may also estimate the 
number densities using equations (I 2 SI 1 and ( ISTT l. 


We apply our diagnostics to the sodium doublet mea- 
sured in HD 189733b using HARPS ground-based data 
dWvttenbach et al.l2015h . At line center, the absorption depths 
and transit radii are 

Fo,i/Fc = 0.9960 ± 0.0007, i?o.i/i?j = 1-227 ± 0.021, 
Fo 2 /F^ = 0.9936 ± 0.0007, i?o 2 /R 3 = 1-277 ± 0.020. 

(35) 

At A A = 0.3 A from line center, the absorption depths and 
transit radii are 

Fi/Fc = 0.9987 ± 0.0007, = 1.168 ± 0.022, 

F 2 IFC = 0.9974 ± 0.0007, = 1-196 ± 0.022. 

All quantities were extracted using Gaussian fits to the data. 
The white-light radius used to normalise the transit radii is 


(1.138 ±0.027)i?j (iTorres, Winn & Holmanl200l . We veri- 
fied that instrumental bro adening is not a major source of error 
dWvttenbach et al.ll2015h . 

To faciliate comparison with iHuitson et al.l (l20T2), we use 
g = 2141 cm s-^ and m = 2.3toh for HD 189733b. We 
propagate the data uncertainities in quadrature. Using the 
isothermal formula, we obtain Ti = 3117 ± 79 K and T 2 — 
4240 ± 103 K. Next, we derive T 1 /T 2 = 0.7284 ± 0.0256, 
which indicates that the non-isothermal treatment is war¬ 
ranted. We obtain 

dT 

— = 0.4376 ±0.0154 Kkm-\ (37) 

oz 

thus confirming the p ositive temperature gradient reported by 
IHuitson et al.l (|2013- Our line-wing and line-center tempera¬ 
tures are 


Ti = 2460 ± 86 K, To,i = 4306 ± 151 K, 

T-u T2,^3336± 117K, ro 2 = 5870A 206K. 

The number densities are 


m = (1.439 ± 0.051) X 10"* cm-3, 
712 = (1-219 ± 0.043) X 10^ cm-3, 
no,i = 3.990 ± 0.140 cm“^, 
no ,2 = 3.128 ± 0.110 cm“^. 


(38) 


(39) 


while b = 13.63 ± 0.48. If the abundance of sodium is solar, 
then the total pressure implied is < 5 nbar. 

The temperatures and temperature gradient inferred using 
our analytical diagn ostics are in broad agre ement with the 
numerical models of iWvttenbach et ^ d2015l). applied to the 
same data set. 


4. CONCLUSION 

We have presented a previously-unelucidated theory for in¬ 
ferring the temperatures, temperature gradient and number 
densities associated with the absorption depths and transit 
radii of the sodium doublet lines. Our formalism is applicable 
to other alkali metal lines such as potassium. Applying our di¬ 
agnostics to infrared lines of molecules requires the inclusion 
of pressure broadening. 
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